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Abstract. Balmer lines are an important diagnostic of stellar atmo- 
spheric structure, since they are formed at a wide range of depths within 
the atmosphere. The different Balmer lines are formed at slightly different 
depths making them useful atmospheric diagnostics. The low sensitivity 
to surface gravity for stars cooler than ~8000 K makes them excellent 
diagnostics in the treatment of atmospheric convection. For hotter stars 
Balmer profiles are sensitive to both effective temperature and surface 
gravity. Provided we know the surface gravity of these stars from some 
other method (e.g. from eclipsing binary systems), we can use them to 
determine effective temperature. 

In previous work, we have found no significant systematic problems 
with using uvhy photometry to determine atmospheric parameters of fun- 
damental (and standard) stars. In fact, uvby was found to be very good 
for obtaining both T^s and log (7. Using Ha and II/3 profiles, we have 
found that both the Canuto & Mazzitteli and standard Kurucz mixing- 
length theory without approximate overshooting are both in agreement to 
within the uncertainties of the fundamental stars. Overshooting models 
were always clearly discrepant. Some evidence was found for significant 
disagreement between all treatments of convection and fundamental val- 
ues around 8000~9000 K, but these results were for fundamental stars 
without fundamental surface gravities. We have used stars with funda- 
mental values of both T^^ and log g to explore this region in more detail. 



1. Introduction 

Balmer lines are an important diagnostic of stellar atmospheric structure since 
they are formed at a wide range of depths within the atmosphere. The depth 
of formation of Ha is higher than that of H/3, thus observations of these profiles 
provide useful diagnostics (e.g. Gardiner 2000). Balmer profiles are relatively 
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insensitive to surface gravity for stars cooler than '--'8000 K (Gray 1992, see also 
Heiter et al. 2002), whilst sensitive to the treatment of atmospheric convection 
(e.g. van't Veer & Megessier 1996, Castelh et al. 1997, Gardiner 2000, Heiter 
et al. 2002). For stars hotter than ~8000 K the profiles are sensitive to both 
effective temperature and surface gravity. However, provided we know surface 
gravity from some other means (e.g. from eclipsing binary systems), we can use 
them to determine effective temperature. 

In previous work, Smalley &: Kupka (1997) found no significant systematic 
problems with uvby and fundamental (and standard) stars. In fact, uvby was 
found to be very good for obtaining Tcff and logg. Using Ha and H/3 profiles, 
Gardiner et al. (1999) found that both the Canuto & MazzitteH (1991, 1992) and 
standard Kurucz (1993) mixing-length theory without overshooting (see Castelli 
et al. 1997) are both in agreement to within the uncertainties of the fundamental 
stars. Overshooting models were always clearly discrepant. However, Gardiner 
et al. (1999) found some evidence for significant disagreement between all treat- 
ments of convection and fundamental values around 8000~9000 K. In this region 
the effects of log g cannot be ignored, and the majority of the T^s stars did not 
have fundamental values of log 51. We have used binary systems with funda- 
mental values of log g, determined fundamental values of T^s and compared the 
results with those from fitting models to Balmer-line profiles. 

2. Effective temperatures of binary systems 

Eclipsing binary systems provide ideal test stars for comparing to models, since 
they enable us to obtain fundamental values of Tcs and log 51. We can obtain 
fundamental values of Tgg provided we know the apparent angular diameter and 
total integrated (bolometric) flux at the Earth. In the case of binary systems, 
where there are no direct measurements of angular diameters, we can obtain 
them from the stellar radius and the parallax of the system (Smalley &; Dworet- 
sky 1995). 

Available spectrophotometry was taken from various sources. Unfortu- 
nately, not all the systems have enough high-resolution spectrophotometry. In 
these cases, however, we have at least UBV or UBVRI colours, which were used 
to estimate the optical fluxes. Near-infrared fluxes were taken from the Gezari et 
al. (1987), 2MASS and DENIS catalogues. The lack of available fluxes for binary 
systems is something that the Citadel ASTRA Spectrophotometer (Adelman et 
al. 2002) could address. 

In most cases the Tegf values obtained using the Hipparcos parallaxes are in 
agreement with that obtained from other methods (e.g. Infrared Flux Method, 
uvhy photometry, etc.) and other determinations (e.g. Jordi et al. 1997, Ribas 
et al. 1998). However, three systems were anomalous and these are discussed in 
detail in Smalley et al. (2002). 

3. Observations 

The Ha and H/? observations were made at the Observatorio del Roque de 

los Muchachos, La Palma using the Richardson-Brealey Spectrograph on the 
1.0m Jacobus Kapteyn Telescope (JKT) in October/November 1997. A 2400 1 
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Table 1. Fundamental values of T^s for binary stars. [Adapted from 
Smalley et al. (2002)] 



Star 






op. 

Tvpcs 


/e 

(10- 


6 Wm^) 


err 






12 Per 


A 


4.16 ± 0.03 


F8V 


195. 


± 12.8 


6371 




176 




B 


4.24 ± 0.03 


G2V 


115. 


± 5.06 


6000 




143 


CD Tau 


A 


4.087 ± 0.010 


F6V 


24.3 ± 3.06 


6110 




415 




B 


4.174 ± 0.012 


F6V 


20.7 ± 1.61 


6260 




397 


UX Men 


A 


4.272 ± 0.009 


F8V 


7.48 ± 1.12 


6171 




302 




B 


4.306 ± 0.009 


F8V 


6.52 ± 0.54 


6130 




233 


13 Aur 


A 


3.930 ± 0.010 


AlV 


2430. ± 250. 


9131 




257 




B 


3.962 ± 0.010 


AlV 


2070. ± 139. 


9015 




182 


WW Aur 


A 


4.187 ± 0.019 


A5m 


58.7 


± 8.65 


7993 




470 




B 


4.143 ± 0.018 


A7m 


49.3 


± 5.01 


7651 




401 


PV Pup 


A 


4.257 ± 0.010 


A8V 


22.0 


± 4.07 


6870 


± 


363 




B 


4.278 ± 0.011 


A8V 


21.0 


± 2.41 


6888 




265 


RS Cha 


A 


4.047 ± 0.023 


A8V 


44.2 


± 5.56 


7525 




307 




B 


3.961 ± 0.021 


A8V 


43.8 


± 3.21 


7178 


± 


225 


HS Hva 


A 


4 326 + 006 


F5V 

J- V 


7.55 


± 0.80 


63Q8 

\JlJiJ\J 




261 




B 


4.354 ± 0.006 


F5V 


6.45 


± 0.40 


6300 


± 


217 


RZ Cha 


A 


3.909 ± 0.009 


F5V 


7.50 


± 1.11 


6440 




444 




B 


3.907 ± 0.010 


F5V 


7.50 


± 0.60 


6440 




396 


7 Vir 


A 


4.21 ± 0.017 


FOV 


978. 


± 80.5 


7143 


± 


451 




B 


4.21 ± 0.017 


FOV 


978. 


± 41.3 


7143 


lb 


433 


DM Vir 


A 


4.108 ± 0.009 


F7V 


3.00 


± 0.53 


5931 


± 


390 




B 


4.106 ± 0.009 


F7V 


3.00 


± 0.28 


5931 




319 


V624 Her 


A 


3.834 ± 0.010 


A3m 


60.0 


± 6.44 


8288 




476 




B 


4.024 ± 0.014 


A7V 


29.0 


± 2.55 


8092 




450 


VI 143 Cyg 


A 


4.323 ± 0.016 


F5V 


56.8 


± 6.32 


6441 




201 




B 


4.324 ± 0.016 


F5V 


53.2 


± 3.38 


6393 




136 


MY Cyg 


A 


4.008 ± 0.021 


FOm 


6.17 


± 0.76 


7459 


± 


891 




B 


4.014 ± 0.021 


FOm 


6.00 


± 0.43 


7408 


lb 


865 


S Equ 


A 


4.34 ± 0.02 


F7V 


210. 


± 16.5 


6393 


lb 


156 




B 


4.34 ± 0.02 


F7V 


210. 


± 8.90 


6393 


lb 


115 
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mm~^ holographic grating was used and a 1124 x 1124 pixel Tek CCD, giving 
a resolution of 0.4A FWHM. Further Ha and H/? observations were made at the 
Mount Stromlo Observatory, Australia in February 2000 using the Cassegrain 
Spectrograph on the ANU 74 inch telescope. A 1200 1 mm~^ blazed grating 
was used, giving a resolution of 0.35A FWHM. Additional H/3 profiles were taken 
from the work of Smalley k. Dworetsky (1995). 

The data reduction of the profiles taken in 1997 and 2000 was performed 
using the Starlink ECHOMOP software package. In most cases the final spectra 
had a signal-to- noise ratio in excess of 100:1. Instrumental sensitivity variations 
were removed from the Ha profiles by comparing to observations of stars with 
intrinsically narrow Balmer profiles, for example early-B or O type stars and G 
type stars, and the H/3 profiles were normalized such that the observed profile 
of Vega agreed to a model with reff=9550 K, log5=3.95, [M/H]=-0.5 (Castelh 
&; Kurucz 1994) and the standard profiles of Peterson (1969). 

4. Effective temperatures from Balmer line profiles 

The observed Balmer line profiles are fitted here to model spectra to compare 
the derived T^s with that from fundamental methods. The following convection 
models were used, using solar-metallicity Kurucz ATLAS models: 

MLT_noOV 1.25 Standard ATLAS9 (Kurucz 1993) models using mixing length 
theory (MLT) without convective overshooting. The value of the MLT 
parameter a is the standard value of 1.25. 

MLT_noOV 0.5 Standard ATLAS9 models using MLT without convective over- 
shooting. The value of the MLT parameter a is 0.5. 

MLT_OV 1.25 Standard ATLAs9 models using MLT with approximate over- 
shooting. The value of the MLT parameter a used is 1.25. 

MLT_OV 0.5 Standard atlas9 models using MLT with approximate convec- 
tive overshooting. The value of the MLT parameter a used is 0.5. 

CM Modified atlas9 models using the Canuto & Mazzitelh (1991,1992) model 
of turbulent convection. 

The synthetic spectra were calculated using uclsyn (Smalley et al. 2001) 
which includes Balmer line profiles calculated using VCS Stark broadening and 
metal absorption lines from the Kurucz & Bell (1995) linelist. This routine 
is based on the balmer routine (Peterson 1969). The synthetic spectra were 
normalized ±100 A to match the observations. The values of Tcfj were obtained 
by fitting model profiles to the observations using the least-square differences. 

Figures 1 &: 2 show the variation of ATeff = Tefj (Balmer) — reff(fund) 
against Tefj(fund) for Ha and H/3, respectively, for the 5 convection models listed 
above. To within the uncertainties, the CM results show no significant variation 
with Teff (fund) for either Ha or H/3. The discrepancy around 8000 K noted by 
Gardiner et al. (1999) is not evident. Even the two anomalous Ha points just 
hotter than 8000 K, for V624 Her, can be brought into agreement if the IRFM 
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Figure 1. Comparison between Tcff calculated from Balmer line pro- 
files Ho; to the Fundamental values. ATeff = Teff (Balmer) - reff(fund) 
is plotted against reff(fund). 
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Figure 2. Comparison between Tcff calculated from Balmer line pro- 
files H/3 to the Fundamental values. ATeff = Teff (Balmer) - reff(fund) 
is plotted against reff(fund). 
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Teff is used (Smalley et al. 2002). The MLT_noOV results are in broad agreement 
with those for CM, but with the a=0.5 models giving better agreement around 
8000 K relative to a=1.25 and CM models. Contrary to Gardiner et al. (1999), 
who reported that F-type stars might require models with a >1.25 (see their 
Fig. 9), we find that the binary systems do not support this. Overall, a=0.5 
models are preferred to those with higher values. The MLT_OV models are 
generally more discrepant, yielding too high values of Tefj (and even larger ones 
for H/3, if a=1.25 is used rather than 0.5), as found previously by Gardiner ct 
al. (1999). Note also the systematic difference between Ha and H/3 for a=1.25 
MLT_noOV models, which is even more pronounced for the MLT_OV models. 



Table 2. A-stars with fundamental values of T^g, but not logff. 



star 


vsim 




CM uvby 


^eff 


Ha 


HP 




km 


Fund 


^eff 


logg 


IRFM 




logs 


TeS logs 


7 Gem 


45 


9220±330 


9250 


3.56 


9040±86 


9220±300 


3.40±0.2 


9060±250 3.52±0.2 


P Leo 


122 


8870±350 


8770 


4.32 


8660±60 


8370±400 


3.77±0.2 


8450±350 4.07±0.2 


a Oph 


240 


7960±330 


7940 


3.80 


7883±63 


7510±100 


3.69±0.3 


7580±150 3.42±0.6 


a Aql 


245 


7990±210 


7840 


4.18 


7588±73 


7420±100 


4.17±0.3 


7450±150 4.38±0.6 


a PsA 


85 


8760±310 


8890 


4.30 


8622±86 


8340±400 


3.87±0.2 





5. The apparent A-star anomaly 

The use of stars with fundamental values of both Tcfj and log 5 has failed to 
support the apparent anomaly around 8000 K found by Gardiner et al. (1999). 
However, there were too few stars within the Tgff range 8000-9000 K to fully 
explore this region. Gardiner et al. (1999) found that four fundamental Tcfj 
stars also showed the anomaly: P Leo, a Oph, a Aql, a PsA. In order to be sure 
that there is no anomaly in the Balmer line profiles, we need to explain why 
these stars might appear anomalous. 

Table 2 summarizes the values of T^s obtained from CM uvby photometry, 
the IRFM and by fitting to Ha and H/3 profiles. We have allowed both T^s 
and log g to vary in order to obtained the best least-squares fit (see Figures 3 & 
4) . Values of log g are also given as obtained from uvby photometry. We have 
also included 7 Gem which is just hotter than 9000 K, but the results are in 
agreement with its fundamental and IRFM T^s values. 

The rapidly rotating star a Aql has recently been studied by van Belle ct al. 
(2001) using interferometry. Their analysis revealed the oblatcncss of the star 
and a new determination of fundamental Teff = 7680=b90 K. This is significantly 
cooler than the previous determination, but in accord with that inferred from the 
IRFM. As such, the Teff from Ha and H/3 are no longer significantly discrepant. 
It is certainly possible that revision to the other fundamental stars could occur 
once new interferometric measurements are obtained, especially a Oph which 
has a similar vsini and might be expected to exhibit significant oblateness. 
Thus, the anomalies for these two stars can be explained in terms of their rapid 
rotation. 

The two other stars, P Leo and a PsA, have lower vsini values, but are 
two most discrepant stars in the Gardiner et al. (1999) sample. Unless the 
fundamental values are truly wrong there must be some other reason for the 
discrepancy. The IRFM values both point to a slightly cooler Teff, but even 
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Figure 3. Ha profiles of A stars. The continuous line is the observed 
profile, while the dotted line is the synthetic profile for best fitting 
parameters given in Table The dash-dot line is that for profiles 
calculated for the fundamental parameters. 
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Figure 4. H/3 profiles of A stars. The continuous line is the observed 
profile, while the dotted line is the synthetic profile for best fitting 
parameters given in Table The dash-dot line is that for profiles 
calculated for the fundamental parameters. 
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then the discrepancy is ^-^300 K. However, in this temperature region the Balmer 
hnes are near their maximum strength and sensitive to logs'. It is certainly 
possible that a small error in adopted logy could lead to a large error in T^s 
obtained from Balmer profiles. In addition, the Balmer profiles change little with 
relatively large changes in Teff. Thus, we conclude that the two stars are not 
discrepant, due to the low sensitivity of Balmer lines with respect to changes in 
Teff and both sensitivity to, and the uncertainty in, the surface gravity for these 
stars. However, it must be noted that the logg obtained from the Balmer lines 
for these stars is systematically lower than that obtained from uvby photometry. 

In general, for stars hotter than 8000 K the sensitivity to logg prevents us 
from using them to obtain values of Tpff to the accuracy required for the present 
task, unless we have accurate fundamental values of log 5. However, until we 
do have stars with accurate fundamental log g values, we cannot be totally sure 
that there is not a problem with the model predictions in this Teff region. 

6. Conclusion 

Balmer line profiles have been fitted to the fundamental binary systems. To 
within the errors of the fundamental T^fi values, neither the Ha or H/3 profiles 
exhibit any significant discrepancies for the CM and MLT without approximate 
overshooting models. As in previous work, the MLT with overshooting models 
are found to be discrepant. Moreover, there are no systematic trends, such 
as offsets, between results from Ha and H/3 as long as a in MLT models is 
chosen small enough (e.g. 0.5). The discrepancies exhibited by the fundamental 
Teff stars in Gardiner et al. (1999) can be explained by rapid rotation in two 
cases and by the fact that the Balmer profiles become sensitive to log g and less 
sensitive to Teg in the other two cases. However, for the time being the lack of 
any stars with fundamental values of both T^s and logg in this region precludes 
the conclusion that there is not a problem with the models in the Tgfj range 8000 
~ 9000 K. 

Full details of this work are given in Smalley et al. (2002). 
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